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The effect of non-magnetic impurities is discussed for both the two-leg Heisenberg ladder sys- 
tem and other spin liquids with excitation gap. It is shown that the random depletion of spins 
introduces a random Berry phase term to the non-linear a model. The classical nature of the 
antiferromagnetic correlation is enhanced by the topological decoherence, and the staggered sus- 
ceptibility shows more singular behavior at low temperatures than the uniform antiferromagnetic 
Heisenberg chain. 
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Recently a class of quantum antiferromagnets which 
has a ground state with a gap in the spin excitation 
spectrum has attracted much attention. These systems 
are characterized either by a local singlet correlation of 
valence-bond-solid state as found in spin-Peierls systems 
[jjj and even-leg spin ladders j^j or by resonating- valence- 
bond (RVB) state as in underdoped cuprates. In such 
systems it was found that the presence of non-magnetic 
impurities replacing magnetic ions often leads to a Curie 
behavior of the uniform susceptibility with a Curie con- 
stant proportional to the impurity concentration. This 
can be understood basically from the idea that, due to 
the depletion of spins, some spins lose their partners to 
form a singlet. They appear, therefore, as nearly in- 
dependent spin doublet degrees of freedom. However, 
the more striking experimental feature is the appearance 
of the antiferromagnetic long-range order (AFLRO) in- 
duced by the non-magnetic impurities. This property 
has been observed essentially in all the systems men- 
tioned above, ||, [l|||] and thus seems to be a universal 
phenomenon in gapped quantum spin liquids. This ef- 
fect has recently been analyzed in terms of the phase 
Hamiltonian for the spin-Peierls system by Fukuyama et 
at U The presence of impurities leads to a local weak- 
ening of the dimerization, which induces enhanced local 
antiferromagnetic correlations. The AFLRO can then 
emerge at low temperatures, because the local antifer- 
romagnetic moments around the impurities correlate in 
phase with each other. The phase Hamiltonian approach 
has also been applied to the two-leg spin ladder system 
with a single non-magnetic impurity. (?j In this case the 
phase field forms a soliton at the impurity position, which 
corresponds to an impurity-induced spin degree of free- 
dom (S = 1/2). Local antiferromagnetic correlation ap- 
pears in the vicinity of the impurity, and would turn 
to an in-phase (quasi) long-range ordered state in the 
ground state. Q Sigrist and Furusaki pointed out that 



low-energy physics of the lightly depleted spin ladder sys- 
tem is described by a spin-^ Heisenberg model with ran- 
dom ferromagnetic and antiferromagnetic couplings. || 
They found that at low temperatures the uniform sus- 
ceptibility follows the Curie law while the specific heat 
shows power-law behavior, C ~ T 2q , where the exponent 
a is expected to be equal to or smaller than 0.25. On 
the other hand, in a recent numerical study Furukawa 
et al. found a significant enhancement of the staggered 
susceptibility x(Q) m t ne presence of non-magnetic im- 
purities. H The asymptotic behavior, however, could not 
be obtained because the size of the system they studied 
is not large enough. 

In this paper we study the effect of non-magnetic im- 
purities on gapped spin liquids via a mapping to the 
non-linear a model with a random topological term. In 
this formulation the appearance of the AFLRO origi- 
nates from the recovery of the topological Berry phase 
term, which suppresses quantum fluctuations such that 
the classical nature of the AFLRO is enhanced. We find 
that the generalized susceptibility x(Q,T) diverges with 
a non-trivial exponent, x(<2,T) ~ j 1 -i-2a^ m ^ ne j ow _ 
temperature limit. This divergence is much stronger 
than the spin-^ antiferromagnetic Heisenberg chain. 
This picture also applies to higher-dimensional systems 
where quantum fluctuations suppress the AFLRO and 
generate the spin gap phase, even though the behavior 
°f x(Q>T) depends on the dimensionality. 

For concreteness we consider the two-leg spin-S* 
Heisenberg ladder as an example of gapped spin liquids. 
We use the functional integral method and write the par- 
tition function Z as 



Z = J DS t (T)e~ 



(1) 



The action A of the system is given by 
A = i5£>({S,(r)}) +JdrJ2 JS i( T ) ■ ( 2 ) 
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The first term is the Berry phase contribution, where 
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u({Si(T)}) is the solid angle enclosed by the closed tra- 
jectory of Si(r) [Si(0) = Si((3)}. @ The second term 
represents the nearest-neighbor antiferromagnetic inter- 
action. 

For the following discussion it is crucial to notice that 
no frustration is induced by the random depletion of 
spins because the bipartite lattice structure remains and 
the antiferromagnetic interaction works only between A- 
and i3-sublattice spins. Consequently we should choose 
the staggered spin as the slowly varying field. On the 
other hand, the Hamiltonian commutes with the total 
spin so that the Fourier components with small wavevec- 
tors are also expected to be slowly varying. From these 
considerations we may express a spin operator in terms 
of two slowly varying fields, fl(xi, t) and 17 (xi, r): 



Si( T ) = sft-iynfa, t) + an( Xi) t) 



(3) 



Here a is the lattice constant, and the two fields satisfy 
the constraint, i? • 17 = and \ fi\ — 1. We can then 
apply the standard procedure to derive the non-linear a 
model. 

Before discussing the randomly depleted Heisenberg 
ladder we first consider pure spin systems. It is well 
known that the effective action A for a single Heisenberg 
chain is ||10| 



1 

2^ 



dr / dx 



I fdf2\ 2 fdf2\ 2 



(4) 

where the dimensionless coupling constant g = 2/S con- 
trols the quantum fluctuations, and the velocity c = 
2JSa sets the energy scale. The Berry phase term ^Bcrry 
is given by 



^Berry = 2tTiQS, 



(5) 



dr 



, n fdO dfi 
ax s 2 • — — x 
V or 



dx 



(G) 



where Q is an integer called Skyrmion number, 
1 ff> 

For half integer S the Berry phase is the multiple of 
7r and thus causes destructive interference between the 
quantum fluctuations with different Skyrmion numbers. 
This leads to the massless spectrum of the triplet excita- 
tions. For integer S, however, the Berry phase plays 
no role for the infinite system and quantum fluctua- 
tions generate the spin gap. In the case of the pure 
two-leg spin ladder, one leg of the ladder contributes 
^■Berry to the effective action, while the other leg gives 
^Berry(^ f2) . The Berry phase terms thus cancel, 
and the action is just the non-linear a model without a 
topological term. |l2[^3) Hence the spin gap opens sim- 
ilarly to the integer spin chain. 

Now we turn to the randomly depleted spin ladder. 
We note that removing spins from the ladder affects the 
effective action in two ways. An obvious effect is the 
local weakening of the coupling J which leads to the 
reduction in the sound velocity c. This tends to decrease 
the energy scale of the problem. The other effect is of 
quantum mechanical nature and is related to the Berry 
phase. The cancelation of the Berry phases as described 



above is violated by the depleted spins, and the Berry 
phase term ^Berry becomes 

^Berry = iS ^ u({S(Xi, t)}) 

= iSj2^({S(xi,r)}) - iS5>({S(x,„,r)}) 

i n 

= -iSj2(-l) in "({n(x in ,r)}), (7) 

n 

where i n specifies the site of a depleted spin. This mod- 
ification leads to a substantial change in the quantum 
mechanical behavior of the system described by eqs. (^) 
and (7), as we will see later. 

In the next step we integrate over the continuum de- 
grees of freedom fi{x ^ Xi n , r) to derive the effective ac- 
tion for Itl{t) = fl{x = Xi n , t). To this end we introduce 
the Lagrange multiplier to implement the constraint 
In(i~) = fl(x = Xi n ,r), and write Z as 



Z= DO DX DIe- Aet! . 



(8) 



The effective action A e s is given by 

AcS = -is]>>irM{/„(T)}) 



c 

"2g 



k 



fl(k,LOl)-fl(-k, -LOi) 



(9) 

where uji is the Matsubara frequency for bosonic fields, 
and (— l) 4 " = +1 (—1) if Xi n is on the A-sublattice (B- 
sublattice). In eq. (9) we have introduced a mass c/lo in 
the second term. This corresponds to taking saddle-point 
approximation for a Lagrange multiplier field represent- 
ing the constraint \fl\ = 1. This can be justified, in the 
limit of small impurity concentration, by noting that l 
is the finite correlation length of the pure two-leg spin 
ladder. We then integrate over fl to obtain 

A eff = -i5^(-l) i "o;({I n (T)}) 

n 

+i y^y^ Xn(-u>i) ■ i n {u>i) 

n uji 

+ 7T ^2 K m,n{ LJ l)^m(^l) ' A n (-W;), (10) 



where we have ignored an additive constant coming from 
integration over i?, and 

mM iJ 7-0O27T W + kl 2^ 

(11) 

In 2 1 1 . For the case where the 
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with h 

concentration 5 of the non-magnetic impurities is very 
small {a/6 3> Iq), we can approximate K^ 1 as 



2k 



'0 ( $n 



i±ie 



(12) 
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After integrating over Xi we obtain, for \u) n \ <C c/Iq, 
A eS = -i5£(-l)*»u,({I n ( T )}) 

71 

~ At V Je-l^-^ +1 '/ io J„(r).I n+1 (r), (13) 

where J = cko/g. Equation (13) can be identified with 
the action for the random-exchange Heisenberg model 
(REHM), 



H — ^ ^ JnSn ' S 



n+1, 



(14) 



where S n = (— l) tn+1 I n - The exchange interaction J„ 
is given by J„ = J(-l) l " +l " +1 exp(-|x l „ - afi^|/Zo), 
which is random in sign as well as in magnitude. || Note 
that the two spins S n and S n +i are coupled ferromag- 
netically (antiferromagnetically) if they are on the same 
(different) sublattice(s). J§| It is important to notice that 
the dependence of the effective action on the impurity 
concentration S appears only in the magnitude of the ex- 
change interaction J n , which is typically of the order of 
Je~ l °/ S . 

The low-energy properties of the REHM have recently 
been studied by Westerberg et al. using a real-space 
renormalization group method. @ [l| They found 
that the correlation length £ for the staggered spin I n 



diverges at low temperatures as 15 



-2<v 



(15) 



where the exponent a was estimated numerically to be 
0.22. Inside a domain of length £ the spins S n align 
almost parallel or antiparallel to each other depending 
on the sign of J n . These strongly correlated spins form 
a single effective spin, whose spin size can be estimated 
from a random-walk argument. Since the sign of 

J n is random, the spin size Stot is given by a random sum 
of positive and negative numbers, thereby yielding Stot ~ 
Consequently, the uniform susceptibility xsCO P er 
site behaves as 



XS(T) 



1 Sf ot 1 
T£(T) T' 



(16) 



and the Curie law persists down to zero temperature. H] 
In contrast to S n , the staggered spins I n align always 
parallel to each other inside the domain. Thus the total 
spin size /tot of the staggered spins I n inside the domain 
is proportional to £. Hence the generalized (staggered) 
susceptibility X/CH behaves as 



Xl(T) 



T£(T) 



rn-l-lot 

T 



(17) 



which is more singular than the staggered susceptibility 
of the spin-i antiferromagnetic Heisenberg chain. This 
implies that the impurity-induced staggered moments J„ 
show at least quasi-long-range order at zero temperature, 
which, even with an infinitesimal inter-ladder coupling, 
will turn into true long-range order, in agreement with 
the experimental observation. 

Let us examine the REHM in eqs. (13) and (Jl4|) in 
more detail. As we discussed above, the spin quan- 



tum number S e s of the effective spin representing the 
correlated spins scales as S e tt ~ ^/£(T), which grows 
with decreasing temperature. This can be considered as 
an intermediate case between the ferromagnetic Heisen- 
berg chain where «S e ff ~ £ and the antiferromagnetic 
Heisenberg chain where Set! is at most S. This means 
that the classical nature of i" is more enhanced in the 
randomly depleted Heisenberg ladder than in the pure 
spin- i antiferromagnetic Heisenberg chain, as we see in 
cq. (p~7]). We summarize these considerations in Table I, 
where four universality classes of one-dimensional quan- 
tum spin chains are classified according to the degree of 
quantum fluctuations of the field I. The Berry phase 
term is most effective in the 1st class while it plays no 
role in the 4th class, where quantum fluctuations dynam- 
ically generate a mass. The even-leg spin ladders without 
impurities also belong to this class. We see that, as the 
Berry phase term becomes more and more effective, the 
classical nature of I emerges, leading to the enhanced 
generalized susceptibility X/(T). It is important to no- 
tice that any small amount of non-magnetic impurities 
is sufficient to change the universality class of the pure 
spin liquids, e.g., the even-leg spin ladders, from the 4th 
class to the 2nd class. This is due to the decoherence by 
the random Berry phase term. A similar idea has been 
discussed recently by Prokof 'ev and Stamp for quantum 
nanomagnets. Jl7|,^8| Here some remarks are in order on 
the random singlet phase, in which spins far apart form 
weakly bound singlet pairs in a random manner. p9|,[20| 
This phase is realized in the random antiferromagnetic 
Heisenberg model where the exchange couplings are all 
antiferromagnetic but are random in their magnitude. 
In this case the Berry phase term is not affected by the 
randomness, and the singlet formation of spins plays an 
essential role. We did not include this phase in Table 
I, because the randomness is not of topological nature. 
Therefore Table I is not intended to exhaust all the uni- 
versality classes of the quantum spin chains. 

Finally, we briefly discuss the two-dimensional case for 
which the mapping to the non-linear a model can be used 
in exactly the same way as in the one-dimensional case. 
For example, in the Heisenberg antiferromagnet on the 
square lattice, Berry phases of each row cancel as in the 
spin ladder, and thus there is no topological term for 
the slowly varying field fi. jl0| For S — \ the dimen- 
sionless coupling constant g ~ 1/5 is not large enough 
for the quantum fluctuations to destroy the AFLRO. 
However, as the mobile holes are doped into the anti- 
ferromagnet (e.g., doped high-Tc cuprates), the AFLRO 
is rapidly suppressed and eventually destroyed. These 
holes are localized when Zn impurities are doped suffi- 
ciently. In this case both the Zn impurities and the local- 
ized holes contribute to the random Berry phase terms 
in the non-linear a model. Then we expect that the an- 
tiferromagnetic correlation is enhanced and the AFLRO 
is recovered with non-magnetic impurities, which is in- 
deed observed in a recent experiment on the underdoped 
cuprates with the (pseudo) spin gap. § The problem of 
the coexistence with the mobile holes and Zn impurities 
is left for future study. 

Before summarizing the results, we point out that a 
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Table I. Universality classes of one-dimensional quantum spin chains classified in terms of the Berry phase. 





class 


Berry phase 


specific heat 


generalized susceptibility 


1. 


fcrromagnets 


iS£ n ^({J n (r)}) 


c ~ Vt 




2. 


random-exchange Heisenberg model 


iSj2 n ^{{IAr)}) (<T„=±1) 


C ~ T 2a 


Xl ^ T -i-2 a 


3. 


antifcrromagnets (half-integer S) 


2niSQ (Q: Skyrmion number) 


C ~ T 


XI ~ T-l 


4. 


antifcrromagncts (integer S) 


2niSQ (Q: Skyrmion number) 


C~e" A / T 


Xj ~ const 



non-linear c model with a random topological term has 
also been discussed in the context of the localization of 
a particle in a random magnetic field. It is an in- 
teresting future problem to explore further this connec- 
tion between these two apparently different problems, 
the random-exchange Heisenberg model and the random 
flux problem. 

In summary, we have studied the effects of non- 
magnetic impurities in gapful spin liquids. We have 
shown that the depletion of spins induces the random 
Berry phases, leading to the destructive interference of 
quantum fluctuations. The effective low-energy model 
for the randomly depleted spin ladder is the random- 
exchange Heisenberg chain, in which the generalized sus- 
ceptibility for the staggered moment diverges as T _1 ~ 2q 
(a w 0.22). The correlation of the staggered spin mo- 
ment becomes more classical due to the decoherence by 
the random Berry phase term. 
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